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Abstract 

We prove the scale invariant Harnack inequality and regularity properties for harmonic 
functions with respect to an isotropic unimodal Levy process with the characteristic ex- 
ponent if: satisfying some scaling condition. We show sharp estimates of the potential 
measure and capacity of balls, and further, under the assumption of that ip satisfies the 
lower scaling condition, sharp estimates of the potential kernel of the underlying process. 
This allow us to establish the Krylov-Safonov type estimate, which is the key ingredient 
in the approach of Bass and Levin, that we follow. 



1 Introduction 



Let X t be a Levy process with the characteristic exponent 

= (x, Ax) - i (x, 7) - / [e l{ - x,z) - 1 - i (x, z) l\ z \<i) v(dz), ze 



^ where A is a symmetric and non-negative definite matrix, v is a Levy measure, i.e. ^({0}) = 0, 
J Rd (1 A \z\ 2 ) v{dz) < 00 and 7 G R d . A generator of this process has the following form, for 

feC 2 b (R d ), 

Af(x) = J2 A ikd%f(x) + J2^d k f(x) + j (f(x + z)-f(x)-l ]z]<1 (z,Vf(x)))v(dz). (1) 

j,k k 

The first exit time of an (open) set D <Z M. d and first hitting time to D c by the process X t 
is denned by the formula 

r D = inf{t > 0; X t <£ D}, T D c = mf{t > : X t E D c }. 
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A function / : R d -»■ [0, oo) is said that is harmonic with respect to X t in an open set D if for 
any bounded open set B such that B C D 

f{x) = E x f(X TB ), xeB. 

The scale invariant Harnack inequality holds for a process X t if for any R > there exists 
a constant C = C(R) such that for any function non-negative on M. d and harmonic in a ball 
5(0, r), r ^ 5, 

sup ^ C inf h(x). 

xeB(0,r/2) ' xeB(0,r/2) 

We say that the global scale invariant Harnack inequality holds if constant in the above inequal- 
ity does not depend on R. 

A measure m(dx) is isotropic unimodal if there exists a non-increasing function m : 
(0, oo) — > [0, oo) such that m(dx) = m (\x\)dx, for x ^ 0. A process is isotropic unimodal 
if a transition probability Pt{dx) is isotropic unimodal, for all t > 0. 

Important class of isotropic unimodal Levy processes are subordinate Brownian motions. 

Let / be a positive function on M. d \ {0}. We say that / satisfies the weak lower scaling 
condition WLSC(/3, 0, C), if /3 > 0, 6 > 0, C > 0, and 

/(Ax) > C\Pf(x), for A > 1, \x\ ^ 6. 

If / satisfies WLSC(/3, 0, C), then we say that / satisfies the global weak lower scaling condition. 

The main purpose of this paper is to prove the scale invariant Harnack inequality and 
regularity properties for harmonic functions with respect to an isotropic unimodal Levy process 
with the characteristic exponent satisfying the weak lower scaling condition. Our main technical 
results are sharp estimates of the potential measure and capacity of balls, and further, under 
the assumption of that if) satisfies the weak lower scaling condition, sharp estimates of the 
potential kernel of the underlying process. This allow us to establish the Krylov-Safonov type 
estimate (see Proposition |2"o"|) . which says that there are c and A < 1 such that for a closed set 
A C 5(0, Ar), 

P X (T A < r B(0 , r) ) > c \B^ r )y X G 5(0 ' Ar) ' 

This estimate is the key ingredient of the proofs of the Harnack inequality and local Holder 
continuity of harmonic functions in the approach of Bass and Levin ([2]) that we follow. 

Our main contribution is the fact that we assume only a mild condition for the characteristic 
exponent. Usually in the existing literature on the Harnack inequality for Levy processes the 
assumptions were given in terms of the behaviour of the Levy measure (see [22], Section 3). 
Our result seems to be important for application to subordinate Brownian motions. There are 
examples when the characteristic exponent is known, while estimates for the Levy measure are 
not. We should also notice that our approach allows to deal with isotropic unimodal processes 
with the Levy-Khinchine exponent behaving at infinity almost like the exponent for a Brownian 
motion, which to our best knowledge were not treated in the literature, except a few particular 
cases. Namely, we can take ip{x) = |x| 2 /(|x|), where I is slowly varying and goes to at infinity. 
An example of such a process is for instance a process with density of its Levy measure equal 
to |x|~ d-2 log~ 2 (2 + |x| _1 ). Moreover, our result allows to extend the scale invariant Harnack 
inequality to its global version for many processes. For instance we get the global scale invariant 
Harnack inequality for a-stable relativistic processes. 

The main results of this paper are following two theorems. The first one is the scale invariant 
Harnack inequality. 
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Theorem 1. Let d ^ 3. If ip satisfies WLSC((3,6,C), then the scale invariant Harnack in- 
equality holds. Moreover, if ip satisfies the global weak lower scaling condition, then the global 
scale invariant Harnack inequality holds. 

The next theorem deals with regularity of harmonic functions 

Theorem 2. Let d ^ 3 and ip satisfy WLSC(/3,9,C). For any R > there exist constants 
c = c(R) and 5 > such that, for any < r ^ R, and any bounded function h, which is 
harmonic in 5(0, r), 

IMz)-%)l<c||/i|| 00 f!^^ , x,yeB{0,r/2). 

Recently there has been a lot of research concerning non-local operators. For instance, 
the paper [BJ established the scale-invariant finite range parabolic Harnack inequality for a 
class of jump-type Markov processes on metric measure spaces. A class of special subordinate 
Brownian motions have been studied in [H], where bounds for the densities of Levy measure 
and potential measure and Harnack inequalities were established. Harnack inequalities and 
regularity estimates for harmonic function with respect to diffusion with jumps are proved in 
[H]. Related work on discontinuous processes include [22], [21], [IS], [I] and [8]. Therefore it is 
pertinent to comment on the differences between our results and those of some related papers. 
For the sake of comparison we present them in the context of Levy processes, however most of 
them are in a more general setting of Markov processes. 

• One of the main assumptions in [6J is that the density of the Levy measure is comparable 
to fnj\\\ x \d on 5(0,1), where / is a strictly increasing continuous function and satisfies 
the fol lowing conditions. There exist < 0i ^ /3 2 , and a constant c such that 

c -i(rA\fM,Jr_A\ 0<ri<r2<oo , 
V r i/ fin) \nj 

2 

s v 

ds < c— — — , r > 0. 



M ^ f(r)' 

One can easily check that the lower scaling condition for / implies the weak lower scaling 
condition for the characteristic exponent, hence our assumption is much weaker than that 
from [6j. In [7], under the assumption that the above estimate for the density of Levy 
measure holds on the whole space, the authors obtained the global parabolic Harnack 
inequality. In our context of isotropic unimodal Levy processes the global lower scaling 
for / is sufficient to get the global Harnack inequality (see Example 2 in subsection 3.4). 

In [22] the following Krylov-Safonov type estimate (Lemma 3.4) was derived 

Such an estimate is sufficient in the proof of the Harnack inequality only if a density 
of Levy measure satisfies similar conditions as in [B]. However, it will not work for 

zy ( x ) = | x |d+2 to%+jipij since aPP^g [t one obtains P X (T A < t b ) > i^j^, for r ^ 1/2. 
Hence if r goes to the term ln vanishes, which makes the above bound useless for 
the proof of the scale invariant Harnack inequality. 
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• In [T4] it was considered a class of special subordinate Brownian motions such that a 
subordinator has a non- increasing density of the Levy measure. Moreover, there was some 
scaling assumption on the Laplace exponent of subordinator in terms of its derivative. 

In the present paper the weak lower scaling condition for the Laplace exponent of the 
subordinator is sufficient to obtain the Harnack inequality and we do not need to assume 
anything else about the Levy measure of the subordinator. This does not mean that 
our result covers all the results of [14J. Their proof is not based on the Krylov-Safonov 
type estimate and it works for a large class of slowly varying Laplace exponents, while 
our approach does not cover that case. This is due to the fact that the Krylov-Safonov 
type estimate does not need to hold for the subordinate Brownian motions driven by 
subordinators with slowly varying Laplace exponents. On the other hand our results 
improve the results from [IB], where it was studied only a particular case of subordinate 
Brownian motion with ip(x) = j^^r^m ~~ 1- 

• Since we do not exclude a case when a Gaussian part is non-zero we mention the paper [8] , 
where diffusions with jumps were considered. In this paper the density of the Levy mea- 
sure is assumed to be bounded from above v(x) ^ c\x\~ d ~ a , for \x\ ^ 1, where a G (0, 2). 
Hence the result can not be applied to the process with A = Id and u(x) = ^ d+2 ln a 1 ( 2 +| a; |- 1 ) • 
Notice that for any Levy process with a non-trivial Gaussian part (rank A = d) the WLSC 
property holds for the characteristic exponent. 

• In (9] it is assumed that, for any r < 1, there exist constants c and a such that v{x — z) ^ 
cr~ a v[y — z), for \x — y\ < r and \x — z\ > r. Therefore for instance this result does not 
cover the case A = Id and u(x) = | a ,| d+2ln ^ 1+ | x |_i- ) e~l a: l , for which we even have the global 

Harnack inequality, due to Theorem [TJ since ip(x) ~ \x\ 2 . 

The paper is organized as follows. In Section 2, we give some preliminary results for general 
Levy processes. Section 3 is devoted to prove estimates of Green function and the main results. 
Moreover, several examples are presented to which our approach applies. In Section 4 some 
conditions are stated that are sufficient to prove the scale invariant Harnack inequality for Levy 
processes not necessarily isotropic and unimodal. 



2 Preliminaries 

In this section we introduce notation and prove some auxiliary results for general Levy processes. 
We denote incomplete Gamma functions by 

t*t /*oo 

7 (<5,*)= / e^u^du, T(5,t)= / e^u^du, 5,t>0. 
Jo Jt 

Let B(x,r) denote a ball of center x and radius r > and let B r = B(0,r). By £ we denote 
the Laplace transform, that means, for a measure /i on [0, oo), 

POO 

Cfi(X) = / e- Xs fi(ds), A ^ 0. 
Jo 

For two non- negative functions / and g we write f(x) ~ g(x) if there is a positive number C 
(i.e. a constant) such that C~ l f(x) ^ g(x) Cf(x). This C is called a comparability constant. 
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We write C = C(a, . . . , z) to emphasize that C depends only on a, . . . , z. An integral J ... we 

understand as f, M 

J [0,6) 

Our primary object is a potential measure G, which is well defined for a transient process, 
by the following formula 

/*oo /*oo 

G(x,A)= P x {X t e A)dt = E x l A (X t )dt, 
Jo Jo 

where A is a Borel subset of M d . In what follows we always consider Borel subsets of R d without 
further mention. Let G(A) = G(0, A). Notice that G(x,A) = G(A — x). By a slight abuse 
of notation we also use G to denote the density of the absolutely continuous (with respect to 
the Lebesque measure) part of the potential measure and then we call G(x,y) = G(y — x) a 
potential kernel. 

The fundamental object of the potential theory is the killed process Xf when exiting the set 
D. It is defined in terms of sample paths up to time tq. More precisely, we have the following 
"change of variables" formula: 

E*f(X t D ) = E*[t<T D ;f(X t )}, t>0. 

The potential measure of the process Xf is called the Green measure and is denoted by 
Gd- That is 

G D (x,A) = E* [ D l A (X t )dt. 
Jo 

The corresponding kernel will be called the Green function of the set D and denoted Gd{x, y). 
If the potential measure is absolutely continuous, then we have 

G D (x, y) = G(y - x) - E x G(X TD - y). (2) 

Another important object in the potential theory of X t is the harmonic measure of the set 
D. It is defined by the formula: 

P D (x, A) = E x [r D < oo;l A (X TD )}. 

The density kernel (with respect to the Lebesgue measure) of the measure Pd(x, A) (if it exists) 
is called the Poisson kernel of the set D. The relationship between the Green function of D 
and the harmonic measure is provided by the Ikeda-Watanabe formula [12] , 

P D (x,A)= [ u(A-y)G D (x,dy), Ac(D) c . (3) 
Jd 

Important examples of isotropic unimodal Levy processes are subordinate Brownian motions 
and some our results are restricted to this class of processes. By T t we denote a subordinator 
i.e. a non-decreasing Levy process starting from 0. The Laplace transform of T t is of the form 

Ee~ XTt = e"^ (A) , A^O, 

where <fi is called the Laplace exponent of T. is a Bernstein function and has the following 
representation: 

0(A) = b\+ [ (l-e- Xu )v(du), 
J (0,00) 
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where b > and \i is a Levy measure on (0, oo) such that J(l A u)\i[du) < oo. 

The potential measure of the subordinator T is denoted by U. Its Laplace transform is 
equal to 

r°° 1 

CU(X) = l e-»U W = m (4) 

We say that a Bernstein function is special if there exists a decreasing positive density u 
on (0, oo) of a measure C/|(o |00 )- F° r a different characterization of special Bernstein functions 
see e.g. [2D| • 

Let _B 4 be a Brownian motion in IR d with the characteristic function of the form 

Ee^ Bt = e~ tlxl \ x G R d . 

By gt(x) we denote the transition density of B t . Assume that B t and T t are stochastically 
independent. Then the process X t = BT t defines a subordinate Brownian motion. It is clear 
that the characteristic function of X t takes the form 

Ee^ Xt = e -^ (|x|2) , x E R d . 

The Levy measure of the process X t is given by the following formula for its density 

POO 

v(x)dx = / g u (x)fi(du)dx, 
Jo 

while its potential measure is equal to 

GiA) = ii^b(A) 1, ° ,( - 4) + jj^w^y- < 5 > 

A subordinator which has a special Laplace exponent <fi is called a special subordinator and the 
corresponding subordinate Brownian motion is called a special subordinate Brownian motion. 

For a function / : -> C we define }'*{u) = sup^^ Re f(x). The following lemma will 
play an important role in the sequel. 

Lemma 3. Let f : M, d — > C be a negative definite function, then 
1 s 2 



2s 2 + 1 



f*(r) </*(sr)< 2(1 + s')f*(r), s,r>0. 



Proof. Since / is negative definite, Re/(x) and / r (x) = f{rx) are negative definite functions 
as well. The upper bound we get e.g. by using [T9] , (1.4) for Re/ r . If s ^ 1, then we get the 
lower bound by monotonicity of /*. For s < 1, by the upper bound 

r(r) = r(rss- 1 )^ 2(1 + s- 2 )r(rs), 

which completes the proof. □ 

Lemma 4. Let f : M d — > [0, oo) and f{u) = sup| x | =n f(x). Suppose that f is positive on (0, oo) 
and f(0) =0. /// satisfies WLSC(/3,6,C), then WLSc(f3,6,C 2 j^ holds for f* , where 

/(Q) _ i 
/•(o) ~ x - 
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Proof. We assume that 9 > 0, since the proof in the case 9 = is similar. Note that, 
WLSC(/3, 9, C) holds for /. Hence, 

f(u) ^ C- X f(a) t 9^u^s. (6) 

Let u ^ 9. Since / is positive on (0, oo), 

/*(«)= max{/*(0), snp /(*)} < snp f( x )^ C - lf -^-f{u). 

0<|cc|<u j (9) e<i\x\^u j (9) 

Hence, for u ^ 9 and A > 1, applying again WLSC(/3, 0, C) for / we arrive at 

f*(Xu) > f(Xu) > C\?f(u) > C 2 j^r\?r(u). 

□ 

To the end of this section we assume that X t is a Levy process characterized by a triplet 
(A, u,i). 

In the proof of following lemma we follow closely the ideas of [23], where authors proved 
similar result for isotropic stable processes. 

Lemma 5. Let D C B r and x G D fl B r / 2 . Then there is a constant C = C(d) such that 

P x (\X TD \ >r) = P D (x,B:) < Ch(r)E x r D , 

where h(r) = \ \A\\r- 2 + \ j + j z [l\ z \ <r - l\z\<\) v{dz)\r~ l + J Rd min(l, \z\ 2 r~ 2 ) v{dz). 

Proof. For / G C 2 (M. d ), by Dynkin formula we have 

G D (Af)(x) = E x f(X TD )-f(x), x G W 1 . (7) 

There is a non-decreasing function g G C 2 ([0, oo)) such that g(s) = 0, for < s < 1/2, 
g(s) = 1, for s > 1. Let c\ = sup s |^"(s)|, then sup s ^ y. We put /(?/) = and 
fr(y) = fiyr" 1 ). Recall that Tr(A) ^ d\\A\\. Hence 

a a2 ft \ ( nt\ w 9'(\y\)\ {y,Ay} g'(\y\) ( d\ 

Since g'(s) = 0, for s ^ 1, 



(*,v/(y)>=^(|y|)^<||z|. 



By for |z| < r, 



2 „ ^ d + 2|d 2 



^r(y, 2) = fr(V + Z) - f r (y) - l\ z \ <r (z, V/ r (j/)) < 4~ SU P d jkfr(v) < C X 

1 v »t 4 r 



2 ' 



j.fc 

And, for |^| ^ r, F r (x,y) ^ 1. By (TjQ) we have 

4/r(y) = / F r (y,z)v(dz) + (j + / (l\ z \ <r -l\ z \ <1 )zu(dz),Vf r (y))+ J^A 3k d%f r (y) 

JR d \ JR d I - k 

< c*±i r ^z) +cl d_±i\m + ^ + ^ (i w <r - i w <x) 



4 J VA<r V" r 2 j 2 r 2 ' 2 
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Applying © to f r (y), we get 

G D (Af r )(x) = E x f r (X TD ), \x\<r/2. (10) 

Since P x (\X TD \ > r) < E x f r {X TD ) and G D 1 = E x t d , the estimates © and (HQJ provide the 
conclusion. □ 

Lemma 6. For any r ^ 7 

8(1 I 2d) ( l|A||r2 + / (lA(r|z|)>(dz)) ^^(r)^2(||A||r 2 + | (l A (r\z\) 2 ) u(dz)^j . 
Proof. Let us observe that 

i/)*(r) ^ sup (z, Az) + sup / (1 — cos (2;, y))v(dy) J ^ 2ip*(r). 

\\z\^r \z\<rJ J 

Since sup| 2 |^ r (z, Az) = \ \A\\r 2 it remains to prove 
773 1 ,v / min(l, (|z|r) 2 ) sup / (1 - cos (z, y))v{dy) ^2 / min(l, {\z\rf) v{dz). 

4(1 + 2dJ J R d \z\^rJ jRd 

Let ^(z) = /(l — cos (z, y))u(dy). Notice that (see e.g. [13], (5.4)), 

(1 - cos((x, y)))g(y)dy, 







2 


1 + 1 


X 


2 



where 



1 



G(y) = - / (2ns)- d / 2 e-^e~ids. 



Hence, by the Fubini-Tonelli theorem 

min(l,(|z|r) 2 ) v{dz) < 2/ - ^L— V {dz 



1 + (|z|r) 2 

2 I / (1 - co8((zr,y)))g(y)dyv(dz) = 2 I ip(yr)g{y)dy. 




Since ip is a negative definite function, by Lemma [3] we have 

/ min(l, (|-2|r) 2 ) v(dz) ^ 4 sup ^(2) / (1 + \y\ 2 )g{y)dy = 4(1 + 2c/) sup $(z). 

JR d \z\^r Jm. d \z\^r 

Since 1 - cosm = 2 sin 2 f < 2 (1 A |m| 2 ), 

/ min(l, |z| 2 |x| 2 ) vidz) ^ --^(x). 

</R d 2 



Hence 



sup ip(z) ^ 2 sup / min(l, (\z\ \x\) 2 ) v[dz) = 2 min(l, (|z|r) 2 ) v(dz), 

\xKr \x\€rJR d JR d 



\x\^r \x\^r . 

which completes the proof of the inequality ffTTT) . □ 
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Since for symmetric processes h(r) = + J ^1 A ^\ vidz), we obtain the following 
corollary. 



Corollary 7. Let X t be symmetric, then 

1 



< /i(r) sC 8(1 + 2d)i)*{r- 1 ). (12) 



Remark 8. Instead of a direct calculation one can compare Pruitt's result [18] and [19] , Remark 
4.8 to obtain ip*^ 1 ) m h(r) under the assumption that there exists a constant c such that 
1 1mip(x)\ ^ cRe^(i), x G M d . 

For subordinate Brownian motions easy calculations improve f[T2"j) . Notice that is increas- 
ing. 

Remark. Let X t be a subordinate Brownian motion, then for r > 0, 

^0(r~ 2 ) < h(r) ^ (1 + 2rf)0(r" 2 ). 

Corollary 9. Let X t be symmetric. There exists a constant C = C(d) such that, for r > 0, 
s < r/2, 

P X (\Xr B(0 J>r)^C^l y \x\<s. (13) 

Moreover if ip satisfies WLSC((3,9,C*), then, there exists a constant C = C (^d,(3,C*,^^ 
such that, for < r < Q~ x , s ^ r/2, 

Proof. Since X t is symmetric, by [19], Remark 4.8, and Lemma El we get 

1 



E x r Ba < E 



where c\ = C\(d). Hence, the first claim follows by Lemma and Corollary while the second 
claim is a consequence of Lemma H] and (|T3|) . We only have to check that if)(u) = s\xp\ x \ =u t/)(x) 
is positive on (0,oo). Suppose that there exists u > 0, such that ^{uq) = 0. Then, by 
subadditivity of y/tf, we have that i/j(nuo) — 0, for any n G N. Hence, by (EJ), ip(x) = 0, for any 
|x| ^ Uq V 9. That implies that ip = 0, what we exclude. □ 



3 Isotropic Unimodal Levy Processes 

In this section we assume that the process X t is isotropic unimodal. In the first subsection we 
obtain estimates for the potential measure and capacity of balls, which are essential for the rest 
of the paper. Next, we use them to get estimates for a potential kernel and Green function 
of the ball. The next subsection contains some improvements of these estimates in the case 
subordinate Brownian motions. Subsection 3.3 is devoted to prove the Harnack inequality and 
regularity estimates for harmonic functions. In the last subsection we give some examples. 
By -00) an d Go we denote radial profiles of ip, v and G, respectively. For instance 

tp(x) = v>o(M). 
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Lemma 10. ([25]) Let X t be a symmetric Levy process, then the following conditions are 
equivalent: 

(1) X t is isotropic unimodal. 

(2) G x (dx) is isotropic unimodal, for A > (in the transient case for A ^ 0), where G x (dx) = 

J °° e- xt P t (dx)dt and P t {dx) = P°(X t e dx). 

(3) A = al, for some a ^ and v is isotropic unimodal. 

Since X t is isotropic its distribution is supported by the whole space, so it is transient for 
d > 3. Notice that G({0}) > iff ^ is bounded. 

In the following proposition we prove that the characteristic exponent of an isotropic uni- 
modal Levy process is almost increasing. 

Proposition 11. We have, for any x G ~R d , 



V>*OD ^ I2ij){x). 



Proof. Let us define, for r ^ 




where Vi(z) = J Rd -i ^(a/M 2 + z2 )dw. Then, we have ipo( r ) = ar 2 + ipo( r )> for some a ^ 0. 
Since v\ is non-increasing on (0, oo), 




We also note that 1 



COSM ^ 7T2U 



2 if \u\ < 7r/3. We have, 





Since the function J Q °° (1 A (zr) 2 ) v\(z)dz is non-decreasing and 1 — cosw ^ 2(1 A u 2 ), 




ip*(r) = sup ip(x) ^ ar 2 + sup 4 / (l A (zs) 2 ) V\(z)dz = ar 2 + 4 / (l A (zr) 2 ) v\(z)d. 



z. 





□ 
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3.1 Green function estimates 

Lemma 12. Assume that d ^ 3. Let f(r) = G [B /f). There exists a constant C\ = C\ 
such that 

C 3fi 

t^v < Cf ^ ^ rFv A > °- 

\?p* (Vx) \?p* 

Proof. Since 



-M 2 1 



( 47T )rf/2 

we have, for A > 0, 



(47r) d / 2 y Rd (4vr) d / 2 7 Rd 

Integrating with respect to dt we have, 



XCf(\) = -\- [ e-N 2 / 4 - ' / ' 



By Proposition HH r= , 0*(|a;|) ^ ip(x) ^ Hence, 

1 . / c -ki 2 /4 ^ < xc f(x) < 12 /" p-i^i 2 / 4 - f/ ' 



(Air) d / 2 J Rd i>*(V\\x\) ^ " (4vr) d / 2 J Rd i>*(V\\x\) 

By Lemma 121 

i f i , f dx 2 y _Mii + r , 

e 4 ^-nrGfe ^ / e 4 _ < _ / e 4 — — - — dx. 



2^*(v / A)^ " 1+M 2 " J Rd V>*(VA>|) ^ ^*(VA) 
The above estimates imply 

C| < A£/(A) < 30 



1 + 1 




2 




;r 


2 



where C x = e ~ " i+jsp^- 

Proposition 13. Assume that d ^ 3. There is a constant C2 = C2{d) such that 



Proof. Let /(r) = G (-B^J. Since f(u) is non-decreasing 

P r°° P 
f(u) ^ - / e- s /«/(^)^ < -C{f){u- 1 ). 

U Ju U 

By Lemma [12] we get 



i/j*(r x ) 

By Lemma |3] we have, for k ^ s, ^,|"_i/ 2 | ^ 4^. Lemma fT2l and ffl~4|) give us, for k > 1, 

11 



where C\ = Hence, for k such that CiT(2, k) — | we have 

< 2 / e- s/u f(s)ds < 2uf{Ku). 
Jo 

Again, by Lemmas [12] and [3] we infer 



2V>*( a /m— 1 ) 4(/t + l)^*(r~ 1 ) ' 

□ 

By Cap A , A ^ 0, we denote the A-capacity with respect to X t . When A = we omit a 
superscript "0". For any non-empty compact set A C M. d there exists a measure p A (see e.g. 
[3]), called the equilibrium measure, such that is supported by A and 



Gp A (F) = J G(F - x)p A (dx) = J P X {T A < oo)dx, F C Mr. (15) 

Moreover p A (A) = Cap(A). If the potential measure is absolutely continuous, then 

Gp A (x) = P X (T A < oo), x G R d . 

Proposition 14. Let d ^ 3. There exists a constant C = C(d) such that, for any r ^ ; 

C-V^r" 1 )^ ^ Cap (Br) < Cty*^- 1 )^. 

Proof. Since d ^ 3, Cap({0}) = 0, so we may assume that r > 0. By Lemma [TU| G is radially 
non-increasing. Let x G £? r , then 



G(y - x)dy. 



G (B r -x)= G({0}) + / 

Jb, 

By [21], Proposition 5.3, there exists a constant cj = Ci(c?) such that 

c a / G(y)dy ^ / G(y - ar)dj/ ^ / G(y)dy. 

J By- J Br J Bf 

Hence 

c x G (K) < G (B~ r - x) <: G (BT) . 

By m, _ _ 

\B r \ = Gp^ (B r ) . 

This preparation yields 

d\B^\ sC G (i^) Cap (B~r) < |^|, 
which combined with Proposition [13] implies 

Cl ' Bl ^*(r-V < Cap fB^) ^ \^A^*(r^)r d . 



36e v ' " v ' y C 2 

□ 
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Remark. Similar calculations provide the following estimates for A-resolvent measure and A- 
capacity 

G\B r ) « A + ^ (r _ 1) and Cap^) « r d (A + r^ 1 )) . 



By [25J, Theorem 1, it is known that the smallest capacity among sets with the same volume 
is attained for a closed ball. Therefore, by the above proposition and Lemma [3] we obtain the 
following result. 

Corollary 15. There exists a constant C3 = C^(d) such that, for any non-empty Borel set A, 

Cap(A) >C 3 r(\A\~ 1/d )\A\. (16) 

To our best knowledge the following upper bound for the potential kernel was known only 
for subordinate Brownian motions and the lower one for special subordinate Brownian motions 
(see e.g. They were obtained as a consequence of appropriate estimates for the potential 

measure and the potential kernel of the subordinator, respectively. 

Theorem 16. Let d ^ 3. Then there exists a constant C4 = C±{d) such that 



G{x) ^ , . A C f. , 1 , , x e R d . 



If additionally ip satisfies WLSC((3, R 1 ,C*), then 

C 

1 1 . , ,, r~ rr < G(x), \x\ ^ bR, 
\x\ d ip*(\x\- 1 ) y ' 11 

where b = c(d, /3)(C*) 1//3 < 1 and C 5 = c(d)b d . 
Proof. Since G is radially non-increasing, 

G{y)dy ^ iB^x^Gix). 



'5(0,1x1) 

By Proposition [13] 



d 36e 
\x\ G(x) ^ 



..Bil^dxl- 1 )' 

which completes the proof of the upper bound. 

Again, by radial monotonicity, we have, for k > 1, 

q, \ > G C^N \ B \x\) > 1 G(B K \ X \\B\ X \) 
" \B K \ X \\B ]X] \ " \x\ d 

Suppose that ip satisfies WLSC(/3, -R -1 , C*), then by Propositions [TBI and fTTl for k\x\ ^ R, 
G (B K \ X \ \ B\ x \) = G (B K \ X \) — G (B\ x \) ^ 



C 2 36e 
ip* ((k|x|) _1 ) 



-1^ 



36e f_ ^(kl -1 ) -A ^ 36e ( C\G\ , 



Hence, for « = (24/(ciC*))?, we get 

36e 1 
v ; |5i|K d |x| d V*(|a;|~ ) 



□ 
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Corollary 17. Let d ^ 3. If WLSC((3, R 1 , C) holds for if), then there exists a constant b < 1 
such that 

G{x)faG{B\ x \)\x\- d , \x\^bR. 

The above comparability is crucial in our proof of the scale invariant Harnack inequality. 
In the next subsection we show the converse of Corollary [T7| in the case of special subordinate 
Brownian motions. 

Remark 18. If if) satisfies WLSC(/3, -R -1 , C*), then a local doubling condition for G holds. 
This means there exists constant C = C(d, (3,C*) such that CG(x) Si G(2x) ^ G(x), for 
< | a; | ^ bR/2, where a constant b is from Theorem [TBI 

Standard arguments provide the following proposition. 

Proposition 19. Let d ^ 3. Suppose that if) satisfies WLSC((3 , R" 1 , C*) , then, for any e G 
(0, 1), there exists a constant L = L(e, d, (3, C*) > 1 such that, for r $C R, 

G Br (x,y) > sG(x - y), L\x - y\ ^ (r - V (r - \y\). 

Proof. Since if) satisfies WLSC, it is unbounded. Therefore G({0}) = and due to Lemma 
[LTj]the potential measure is absolutely continuous. We may and do assume that \y\ ^ |x| and 
L ^ b^ 1 . Since \X TBr — y\ ^ r — \y\ ^ L\x — y\, by radial monotonicity of G and (j2J), 

G Bp (a;,y) = G(x - y) - E x G(X TBr - y) > G(x - y) - G(L(x - y)). 

By this, Theorem [TS] and Lemma [3J 

™ / x ™ , . / C 4 ^*(lx-d- x ) \ ^ s ( . 4C 4 \ 

G Br (x,y) > G(« - ,) - ^ / (( ^ _y 1} j > C(« - 1/) (l - ■ 

Hence, for L = ( Cr ^l e ) J V 6 we obtain 

G Br (x,y) > zG(x - y). 

□ 



3.2 Subordinate Brownian motions 

In this subsection we improve Theorem [161 and Proposition [19] in the case of subordinate Brow- 
nian motions. Namely, we prove that 6 = 1 and L = 2, for some e > 0. We assume in this 
subsection that X t is a subordinate Brownian motion. 

The following lemma is well known (see e.g. [ID]), but for the convenience of the reader we 
prove it with a short and simple proof. 

Proposition 20. For r > 

— £ U\0,r) ^ 
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Proof. Notice that for A > 1, <ft(\r) ^ A0(r). Hence, 



20((2r)- 1 ) 0(r- 

/•oo 

^ U[0,2r) + e- l J^ e-^~ H U{dt) < C/[0, 2r) + ( h 



e- 1 



which proves the lower bound. 
On the other hand 



U[0,r) < e / e~ r ~ H U{dt) < ' 



r -iy 



□ 



The following theorem is an improvement of Theorem [161 Such result is known (see e.g. 
|26j . Theorem 1), under an additional assumption that is a special Bernstein function. 

Theorem 21. Let d ^ 3 andX t be a subordinate Brownian motion. Iftfi satisfies WLSC(j3, R~ 2 , C*), 
then there exists a constant C 6 = C 6 (d, (3,C*) such that 

Proof. Let k < 1. By (JSJ we have 

,1*13 

^ / g t {x)U{dt) > {g K {l) A ^(1)) ^-^laf, |x| 2 ), 

A|:r| 2 

where 1 = (1,0, ...,0). Suppose that satisfies WLSC(/3, i?" 2 , C*), then by Lemma l20| for 
|z| ^ i?, 

f/[/t|x| 2 , |x| 2 ) = t/[0,|x| 2 ) -f/[0,/t|x| 2 ) ^ ~ 



20(M~~ 2 ) 0((k|x| 



1 - 2e " - C1 ^ > 1 - 2e " (1 - CV) 



2(f)(\x\- 2 )\ 1 ^((nlxl 2 )- 1 ) J 2<f)(\x 



where c 2 = ( e _2)c* • Hence, for k = (2c 2 ) , we get 



|x| rf 0(|x| 



wherec 3 = ^(^(l)A^(l)). □ 

The following theorem is a converse of the above theorem (and Corollary [TT]) in the case 
of special subordinate Brownian motions. Since the comparability in Corollary [T7] is the key 
ingredient in the proof of the Krylov-Safonov estimate it seems that the approach of Bass and 
Levin for proving the Harnack inequality can not be used if does not satisfy WLSC. 



Theorem 22. Let d ^ 3 and X t be a special subordinate Brownian motion. There exists a 

c 

\x\ d <f>{\x\- 2 )- 



constant C such that G(x) ^ \x\dM\x\- 2 ) > f or N ^ R $ ' $ sa ^ s fi es WLSC((3, R 2 , 1), for some 



(3>0. 
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Proof. Due to Theorem [5T] it is enough to show that the existence of a constant C\ such that 

Cl 



G(x) > 



\x\ d (j){\x\ 



(17) 



implies the weak lower scaling condition for 0. Suppose that (JTTj) holds. Since the process is 
transient 



I G(x)dx ^ G{B R ) < oo, 
Jb r 



which combined with (1171) shows that <f> is unbounded and consequently the potential measure 
of X t is absolutely continuous. Since is a special Bernstein function, by (jHJ), we have 



G{x) 



g s {x)u(s)ds, 



where u is non-increasing. By 
1 



(A) 



se- Xs u{s)ds ^ m(A _1 ) / se- Xs ds = (1 - 2 e - 1 )A~y A" 1 ). (18) 



Since i is completely monotone 
and w implies 



is non-increasing. Due to (jl8p . monotonicity of 
1 



1 - 2e 



-i 



v s - 2 ) . 



Hence, 



lV 



I — 2> 



1 - 2e-! 

where c 2 = ^g^M^il . Using 03 



(|x| 4 V s 2 ) g s (x)ds = c 2 



-i)'(\-n 



\x 



d+2 



C 2 



j(A) < ^ (-^) (A)A, for A ^ /T 2 . 

That is, for A ^ -R -2 , 0(A) ^ ^-(//(A)A. Let < ^, then the function (f){u)u~^ is increasing on 
[i?~ 2 ,oo), since (0(M)n _/3 )' > 0. In consequence 

0(Am)(Am)-^ 



> 1, 



A > l,u > iT 



which completes the proof. 



□ 



Let D be an open set and x E D. Denote Su(x) a distance x from a boundary of D. The 
following theorem improves Proposition [JJB Like in the case of Theorem [5T] such result was 
known only for special subordinate Brownian motions for which the characteristic exponent or 
its derivative satisfies some scaling conditions (see e.g. [IS], [II])- These results were obtained 
by standard arguments we used in Proposition [T9l therefore the appropriate constants depend 
on a process. Our proof for a special subordinate Brownian motion does not require any scaling 
properties and the appearing constant depends only on the dimension. 
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Theorem 23. Let d ^ 3 and D be an open set. Suppose that <fi is a unbounded special Bernstein 
function, then 

G D (x,y)^C 7 G(x-y), 2\x - y\ ^ 5 D (x) V 8 D {y). (19) 

where C 7 = T -^0± (l - e-'iy 

If <fi is only a Bernstein function but satisfies WLSC((3, R~ 2 ,C*) , then ( Q7| ) holds, if addi- 
tionally \x — y\ ^ R with a constant C — C(d, (3, C*) instead of C 7 . 

Proof. Let us assume that 5d{x) ^ 5d(v) and x ^ y. Since 2\x — y\ ^ b~v{y) ^ \X TD — y\, by 
(J2J) and radial monotonicity, 

G D (x, y) = G(x, y) - E X G(X TD - y) > G(x - y) - G{2{x - y)). (20) 

Let us define a function 

K p (q) = (47Ts)- d/2 (e-£ - e'ij U(ds), p,q^0. 
Due to §5§ we have 

G(x -y)- G(2(x - y)) = K\ x „ y \{2\x - y\). (21) 

Moreover, 



f\ X ~y\ 2 , l-._-.l2 / 3\ x -v\ 2 \ 

K\ x . v \(2\x-y\) > J (4vr S )- d /V^ M - e"^ J C/(ds) 

/ , x r\x-v\ 2 , |__„,2 

> (l-e-^J y (47rs)- d/2 e-^f/(rfs). (22) 

Suppose that is a special Bernstein function. Then there exists a non-increasing function 
u such that U(ds) = u(s)ds. Monotonicity of u yields 

2 n_| n 

<*/2 -^^,f W <T J__ /, _„|2N|_ _„|2-_ 



(47r S )-^e-^ M ( S )^ < ±1-—LAL U (1 X - y \*)\ x - y 



, j 47T d / 2 



r (| - 1, J) 

Hence, 



7 (_ _ 1 I) f k-y| 2 , . _ ,2 

< ^ /t / {^s)-^e-^u{s)ds. 

1 U _ ^7 JO 



G(x - y) < ^1 + r(!_i I) ) jf ^ ( 47rs )^ /2e "^ M ( s ) rfs < Pa^l(2|* - 2/1), 
r(f-i) 

where c 2 = — x , — jv. In consequence 

r(|-i,*)(i-«-*) 

G D (x,y)> c^ l G{x -y), 2\x - y\ ^ 5 D (y). 

Now, let us suppose that is only Bernstein function (it is no longer assumed that it is 
special) satisfying WLSC(/3, R~ 2 , C*). Then, by the proof of Theorem I2T1 and Theorem [161 there 
exists a constant k < 1, such that 

\x-y\ 2 

g s (x - y)U{ds) > — G(x -y), \x - y\ < R. 

K\x~y\ 2 W 
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Hence, by (j2Dj)-(j22j) 

G D (x,y)^ (\- e -l)Q-G{x-y). 

□ 

Remark. Let d ^ 3. Suppose that <ft is an unbounded special Bernstein function then there 
exists a constant C = C(d) such that, for any r > 

CG(x -y)^ G Br (x, y) ^G(x-y), x, y e B(0, r/5). 

If cf) is only a Bernstein function satisfying WLSC(/3, -R -2 , C*), then there exists a constant 
C — C(d, (3, C*) such that, for any r ^ R the above inequality holds. 

3.3 Harnack inequality and Holder regularity 

The goal of this subsection is to prove the main results of this paper, that is Theorems [T] 
and [2J In this subsection we assume that X t is an isotropic unimodal Levy process with the 
characteristic exponent satisfying WLSC(/3, 9, C*). Since ip satisfies the weak lower scaling 
condition, therefore is unbounded. Hence, the potential measure is absolutely continuous. Let 
R = O^ 1 . By L we denote the constant from Proposition [19] for e = Cj or L = 2 in the case of 
special subordinate Brownian motions and let r = jLTI- 

In the proof of the following proposition we follow closely the ideas of jl], where symmetric 
stable Levy processes were considered. 

Proposition 24. Let d ^ 3 and ip satisfy WLSC((3, R" 1 ,C*). Then there exists a constant 
Cg = C(d, /3,C*) such that, for any r ^ R , and any non-negative function H such that 
supp H C B r . , 

E X H(X TB J ^ C % EPH(X TBr ), x,y e B r 

Proof. Due to Lemma [TO] the Ikeda- Watanabe formula ([3]) and (T2J) we obtain that the Poisson 
kernel of B r exists and 



E»H(X TBr ) = f H(z)P B , 

JBZ 



(y,z)dz, 



where 



P Br (y,z)= G Br {y,w)u(z -w)dw. 

J B r 

Hence, it is enough to prove there is a constant c\ such that 

PB ro (x,z) < dPsXViZ), x,y E Bzo., \z\ > r. (23) 
By Proposition [19] and radial monotonicity of G, 

P Br (y,z)> [ G Br (y,w)v(z-w)dw>C 7 G (2r ) j v{z-w)dw. (24) 

J Brr. J B rt , 



Since v is radially non-increasing, for w £ B\ 



ro ' 



Ib( w 20 v ^ u ~~ z )d u 

v[w - z) < — {Wz ' 8 / < c 2 r Q d I v[u - z)du, 

\B r 1 



'12. J R 

» I J -Org 
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where w z = w + anc ^ c 2 = ]Jm- Hence, by a doubling condition (see Remark ITS]) and 

radial monotonicity of G, 



P Br (x,z) ^ / G{x — w)u(w — z)dw 

JB ro 

^ G ( -j J / f(w — z)dit; + c 2 r^ d / z/(w — / G(x — w)dw 

^ C3Go(2ro) + C2rQ d / G(x — w)dw J / v{u — z)du 

< c 4 (G (2r ) + (2r )- d G(B 2ro )) f u(u - z)du, 

J B r(j 

where c 4 = c 3 V (2 d c 2 ). Corollary [T7I provides 

Pe ro (x, z) ^ c 5 G (2r ) / i/(iz - z)du, 

J B T q 

for some constant c 5 = c 5 (<i, (3, C*). Due to (J2"l|) this implies f l2"3"j) . which completes the proof. □ 

In the following proposition we prove the Krylov-Safonov estimate, which is the crucial for 
proving the scale invariant Harnack inequality. In the proof we use some ideas of [21], Lemma 
6.2. 

Proposition 25. Let d ^ 3 and ip satisfy WLSC((3, R~ l ,C*) . There exists a constant Cg = 
C(d, (3, C*) such that for any r ^ R and any compact A C B ro , 

P x (T A <r Br )^C 9 ^-, xeB ro . 

Proof. Let B be an open set and A C B be compact. Similarly to (II 5p . let 

GbPa(x) = / G B (x,y)p A (dy). 



A 



Then, by the strong Markov property and (J2J), 

GbPa(x) = G PA (x) - E x Gp A {X TB ) = P X (T A < oo) - (T A < oo) 

*C P X (T A < oo) - E x [P x tb(T a < oo),T A > r s ] = P X (T A < r B ). 



On the other hand 



This implies 



G B p A (x) ^ MG B (x,y)p A (A). 

y&A 



P X (T A < t b ) > inf G B (x, z) Cap(A). (25) 



By Proposition [19] and radial monotonicity of G, for s £ i? 



inf G Br (x,z) ^ C 7 inf G(x - z) ^ CVGo(2r ) 

z£B rQ z£B rQ 
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By Theorem Q2J G (2r ) ^ (2ro) d r % ro) -i) • Combining this with (l2"5"j) for A C B = B r , and ([To]) 
we obtain 

where c\ = CiC ^f' 1 . By Lemma El for A C B ro , there exists a constant c 2 = c 2 (cf) such that 

V>*((2r )- 1 ) * °2- ■ nt,11Lt; ' 



P x (T A <T Br )^ Cl c 2 \B 



□ 

Let us notice that until now, under the assumption that X t is isotropic unimodal with its 
characteristic exponent satisfying WLSC(/3, -R -1 , C*), all the constants that appear in the paper 
depend only on d, (3, C*. None of them depends on R or 9, respectively 

Now, we are ready to prove the main results of our paper. 

Proof of Theorem [IJ We prove the result for bounded harmonic functions. The boundedness 
assumption can be removed in a similar way as in [22], Theorem 2.4. Assume that ip satisfy 
WLSC(/3, 9, C*). Let R > 0. We prove that there exists a constant c\ = ci(R ) such that, for 
any function h non- negative on M d and harmonic in a ball B r , r ^ R , 

sup h{x) sC C inf h{x). (26) 

xeB(0,r/2) ' x£B(0,r/2) 

Recall that R = Q~ x . With Propositions I!M1 and [231 at hand we can use the approach of Bass 
and Levin (|2J) to get the existence of constants c 2 = c 2 (d, (3,C*) and a = a(d,(3,C*) < 1 such 
that, for any function h non-negative and bounded on M d and harmonic in a ball B r , r ^ R, 

sup h(x) ^ c 2 inf h{x). 

X£B ar x£B ar 

Next, we use the standard chain argument to get 

sup h(x) ^ C3 inf h(x), 

x&B r/2 x£B r/2 

where c 3 = c 3 (d, c 2 ,a). If Ro ^ R we have fl26|) . Notice, that if t/> satisfies the global weak 
lower scaling condition (R = 00) we get the global scale invariant Harnack inequality, since 
we can take C\ = c 3 and c 3 does not depend on R . For R > R, one can use again the 
chain argument to get (12"6"1) . for any harmonic function on B r , r ^ R . But then the constant 
Ci = ci (c/,c 3 ,f ). □ 

To deal with dimension d < 2 we use the idea from [T7] , which relies on extending harmonic 
functions to higher dimensional spaces. 

Corollary 26. Let d ^ 2. Suppose that there exists an unimodal isotropic Levy process Y t G M 3 , 
such that X t is a projection ofY t . Ifip satisfy WLSC(f3, 9, C*), then the scale invariant Harnack 
inequality holds. 

Proof. We present only the one- dimensional case. Without loss of generality we can assume 
that X t = Y^, where Y t = iY^ l \ Y^ 2 \ Y^). Suppose that h is harmonic and non-negative with 
respect to X t in (— r, r), then by the strong Markov property a function / : IR 3 — > [0, 00) defined 
by f(x^\ x*- 2 -*, x^) = h(x^) is harmonic with respect to Y t in (— r, r) x M 2 . Since X t is isotropic 
we obtain that the characteristic exponent of the process Y t satisfies WLSC(/3, -R -1 , C*). Due 
to Theorem [T] the scale invariant Harnack inequality holds for Y t , so it must hold for X t . □ 
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Proof of Theorem [H With the Krylov-Safonov type estimate (Proposition |2"51) and the second 
part of Corollary [9] the proof is similar to the proof in [2], Theorem 4.1, therefore it is omitted. 

□ 

Corollary 27. Let d ^ 2. Suppose that there exists an unimodal isotropic Levy process Y t G M 3 , 
such that X t is a projection ofY t . Ifip satisfy WLSC(/3, 8, C*), then the conclusion of Theorem 
d holds for X t . 

Let us remark that in general we can not find an isotropic Levy process Y t in higher di- 
mension such that X t is a projection of Y t . If there exists such process Y t , then ^od^l) = 
J K 3-d v o {\/\ x \ 2 + \y\ 2 )dy- Hence z/ must be continuous on (0,oo). Therefore if is is not con- 
tinuous the construction of Y t is impossible. 

On the other hand any Bernstein function defines a subordinate Brownian motion in every 
dimension hence the following theorem holds with no restriction on dimension. 

Theorem 28. Let d ^ 1 and X t be a subordinate Brownian motion. Suppose that satisfies 
WLSC(j3,6,C*). Then the scale invariant Harnack inequality as well as the conclusion of 
Theorem^ hold. Moreover, if satisfies the global weak lower scaling condition, then the global 
scale invariant Harnack inequality holds. 



3.4 Examples 

Recall that u(dx) = v$(\x\)dx. If we do not mention d ^ 3. 

Example 1. Let A = and vo(r) = r G (0,1), where f{r) is non-increasing and non- 
negative. If /(Ar) ^ c\~Pf(r), for A > 1 and Ar ^ 1 then the scale invariant Harnack inequality 
holds, since ip satisfies WLSC(/3, 9, C) for some ^ and C > 0. 

Example 2. Let A = and vo( r ) — ^r, r £ (0, oo), where /(r) is non-increasing and 
non-negative. If /(Ar) ^ cA _/3 /(r), for r > and A > 1 then the global scale invariant 
Harnack inequality holds. For instance this example is applicable for the following processes, 
(«,«! G (0,2)): 

• Relativistic stable process (f(r) ~ r _a (l + r Y a+d ~ 1 ^ 2 e~ r ), for d ^ 1. 

• Truncated stable process (f(r) = r _a l(o i i)(r)). 

• Tempered stable process (f(r) = r~ a e~ r ). 

• Isotropic Lamperti stable process (f(r) = re Sr (e s — l) - "™ 1 , 5 < a + 1). 

• Layered stable process (f(r) = r~ a l( 0i i)(r) + r~ ai l[ lj00 )(r)). 

The scale invariant Harnack inequality for all these examples are known, for instance by [6], 
but to our best knowledge the global one only for the last one (see e.g. [7]). Another example 
to which our result applies is f(r) = r~ 2 log _2 (l + r -5 ), for 5 < 1. Note that / does not satisfy 
the condition (1.5) in [6], so the scale invariant Harnack inequality can not be concluded from 

m- 

Example 3. Let be a Bernstein function comparable with a function regularly varying 
at infinity with index a. If a G (0, 1] the scale invariant Harnack inequality holds for the 
corresponding subordinate Brownian motion. This covers for instance results of [16], where a 
particular 0(A) = log (x+A) — ^ was cons idered, for which we even have the global scale invariant 
Harnack inequality due to Theorem [TJ 
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Example 4. Let satisfy WLSC(/3i, 0, c x ) and Vi satisfy WLSC(/3 2 , 0, c 2 ). Then Vi + V>2 
satisfies WLSC(/?i A /3 2 , 0, ci A c 2 ). Hence, if ipi, 1P2 are Levy-Khinchine exponents of isotropic 
unimodal Levy processes, then the global scale invariant Harnack inequality holds with constant 
depending only on dimension, (3\ A /3 2 and c\ A c 2 . In particular the global scale invariant 
Harnack inequality holds for a sum of two independent isotropic a-stable process with exponents 
if>i(x) = bi\x\ ai and ?/> 2 (x) = 6 2 |x| Q2 , where < a.\ ^ a 2 ^ 2. Moreover, the constant in the 
Harnack inequality depends only on dimension and a\ in this case. 

Example 5. Let X t be an isotropic unimodal Levy process with the characteristic exponent 
if), independent of a Brownian motion B t , then the scale invariant Harnack inequality holds for 
X t + aB t , a > 0. If additionally if) satisfies the following 

if)(x) ^ Ck~^iP{kx\ \x\ ^ 1, k < 1, 

for some constants C and (3 > 0, then the global scale invariant Harnack inequality holds for 
X t + aB t , a > 0. 

Of course for all of the above examples Holder continuity for bounded harmonic functions 
holds as well. 



4 Applications to more general Levy Processes 

Let X t be a general Levy process and d ^ 3. In this section we relax the assumptions and 
comment on validity of the previous results in this new setting. 

We set three conditions which to some extent replace the core assumption of the previous 
section that the process is isotropic unimodal. 

(Al) Assume that v(dx) = v{x)dx and there exist constants C*, R > such that 



v[x — y) ^ C*r j v{y — z)dz, for any r <\x — y\/2 A R. 

JB(x,r) 

(A2) Assume that G(dx) = G(x)dx, x 7^ 0, and there are constants C\, R > such that 
(C 2 ) _1 G(|x|) ^ G(x) ^ C 2 G(|x|), for |x| ^ R and G is non-increasing. 

(A3) There exists a constant C3 such that | lmif)(x)\ ^ Cg Reip(x) and ip*(\x\) ^ Q Reif)(x), 
x E R d . 

Notice that under (A3) process is transient (d ^ 3). 

In Remark [8] we explain that the claim of Corollary [7] holds if | Im -0 (a;) | ^ lie if) (x), 
x G R d . Of course then the comparability constant will be depend on C3. This condition is 
also sufficient to get f[T3"j) . The second claim of Corollary M holds if we assume additionally 
that Re?/' satisfies WLSC(/3, 8, C). If we assume (A3) we infer the claim of Lemma [T2"l Indeed, 
under (A3) we have 

\Cf(x)- 1 I c -H 2 /4 Re^(x) ^ C |,| 2/4 dx 

In the proof of Proposition [13] we used only Lemma [12] and Lemma [3] hence the conclusion of 
Proposition [13] holds under (A3). 

In the proof of a counterpart of Proposition and f lTB]) we use the following theorem. 
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Theorem 29. ( 7??]/ . Theorem 3.3) Let X\ and X 2 be Levy processes having exponents ip\ and 
ip2, and capacities Cap* and Cap* respectively. If X > and C > are such that 



Re — — < C Re — — for all x, 

then 

Cap* (A) < AC Cap* (A) 

for any analytic set A. 

Lemma 30. For any Levy process X t there exists a subordinate Brownian motion with the 
characteristic exponent <p(\x\ 2 ) such that, for r ^ 0, 

(r 2 ) ^ 4j*(r) < 40(r 2 ). 



8(1 + 2d) 

Proof. Let us define a Bernstein function by the formula 



/ (l -e~ |z|V ) v{dz) + \\A\\r. 



Since § (1 A u) «C 1 - e~ u ^ 1 A u 



(1 A (|^|r) 2 ) i/(cfe) + ||A||r 2 ^ 0(r 2 ) ^ / (l A (|^|r) 2 ) + \\A\\r 



1 
2 

which completes the proof due to Lemma □ 

The following proposition is a counterpart of Proposition \\M and (TL6]) . 
Proposition 31. Let (A3) hold. Then 

Cap(B r ) ^ rV(r ! ). 
Moreover, for any non-empty Borel set A we get 

Cap(A) > C(d,C;)r(\A\- 1/d )\A\. 
Proof. Let A > 0. Then, by (A3), 

SRe f_i\< c- 3 



(1 + (Q) 2 )(A + V*(M)) \A + ^(x) / A + ^*(|ac|) 

Hence, by Theorem I2U1 and Lemma I3T)1 we obtain, for any analytic set A, 

32(1 + 8d)C* 3 CaP ^ (A) ^ CapA(A) ^ 16(1 + {C;)2) Cap ^ (A) ' 

where Cap^ denote the capacity of a subordinate Brownian motion Y t with the characteristic 
exponent <fi defined in Lemma 13171 Since Cap (A) = lim A _ s>0 + Cap x (A) the above inequality holds 
also for A = 0. Finally, we use (TH|) . ffTEj) for F t and again Lemma 13111 to get the conclusion. □ 
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To get conclusions of Theorem (THl Corollary [T7] and Proposition [TH] it is enough to assume 
(A2), (A3) and the weak lower scaling condition for Re0. Under the same assumptions Propo- 
sition |25] holds. We additionally need to assume (Al) to prove Proposition [2H Finally we have 
the following theorems. 

Theorem 32. Suppose that (A1)-(A3) hold and Reip satisfies the weak lower scaling condition, 
then the scale invariant Harnack inequality holds. 

Theorem 33. Suppose that (A1)-(A3) hold and Reip satisfies the weak lower scaling condition, 
then the conclusion of Theorem [H hold. 

Remark 34. If (Al), (A2) holds for R = oo and Reip satisfies the global weak lower scaling 
condition then the global scale invariant Harnack inequality holds. 

For instance, we can use our results to the sum of two independent isotropic stable processes 
with drift. More precisely we consider a process X t with ip(x) = \x\ ai + \x\ a2 — i (2,7), where 
1 < ct\ < 2, a 2 ^ 1 and 7 G M. d . It is easy to see that this process satisfy (Al) and (A3) and 
the global weak lower scaling condition. Since estimates for the heat kernel of this process are 
known locally in time (see |24j) and estimates for the heat kernel of the sum of two independent 
isotropic stable process ([?]) one can check that potential kernels of these two processes are 
locally comparable. Hence the assumption (A2) is satisfied. Therefore we infer that the scale 
invariant Harnack inequality holds for X t . 

Note that similar condition as (Al) appeared in [I] and [Tj and exactly the same in [5] and 
[5J. Instead of conditions (A2) and (A3) the authors of the above mentioned papers assumed 
some additional conditions for a Levy measure. 
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